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Abstract 

We propose a new method to investigate interactions involved in atomic force 
microscopy (AFM). It is a dynamical method relying on the growth of oscil- 
lations via parametric resonance. With this method the second and third 
derivatives of the tip-surface interaction potential can be measured simulta- 
neously. Because of its threshold behavior parametric resonance AFM leads 

to sharp contrasts in surface imaging. 
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Atomic force microscopy (AFM) is a local method (on the scale from nanometers to 
microns) able to sense the interaction potential between a relatively fiat surface and a more 
or less well-defined probing tip mounted on a cantilever [|I|. 

The different methods to sense the interaction potential can be divided into two groups: 
Static methods where the lever deflection proportional to the force between tip and surface 
is measured and dynamic methods where the frequency shift of the harmonically driven 
lever is measured p|. From static methods one can thus get the first derivative of the 
interaction potential as a function of the tip-surface distance. With dynamic methods where 
the oscillation amplitude is held small one gets the force gradient, i.e., the second derivative of 
the interaction potential by measuring the shift of the resonance frequency. There are also 
dynamical methods with large oscillation amplitudes which involve nonlinear oscillations 
induced by the interaction In all these dynamical methods the AFM cantilever is 

excited near of one of its resonance frequencies. 

In this Letter we propose a new dynamical method which is able to measure the second 
and third derivative of the interaction potential simultaneously. Contrary to the dynamical 
methods described above, the AFM cantilever is excited near a frequency which is roughly 
twice its resonance frequency. First we discuss qualitatively the main properties of such 
a parametric resonance. In particular, we will see from which kind of measurable data 
one can obtain the third derivative of the potential. Following this qualitative discussion a 
single-oscillator model is treated quantitatively. This leads to convenient (but approximate) 
formulas for the second and third derivative. Finally, we propose an operational scheme for 
measuring the quantities from which one obtains the derivatives in question. 

The reader might ask, can anything interesting happen by driving an AFM cantilever 
away from its resonance frequency? The answer is: It is possible to rapidly detect the growth 
of an instability due to the nonlinear tip-surface interaction. The instability mechanism is 
the first-order parametric resonance 0. Parametric resonance occurs in an oscillator which 
is driven by a modulation of its eigenfrequency rather then by an external force. The simplest 
example is a pendulum with an oscillating length. Such oscillations can be excited very easily 
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if the modulation frequency is twice the average eigen frequency of the oscillator. Everybody 
knows this intuitively from the experience of a swing-boat at a fair. 

How does this mechanism manifest itself in an AFM cantilever which is driven in the 
usual way? Upon approaching the sample, the resonance frequency of the cantilever changes 
slightly by an amount proportional to the second derivative of the interaction potential. 
Thus the resonance frequency of the cantilever depends on the tip-surface distance. If the 
cantilever is driven out of resonance, it nevertheless oscillates with an amplitude which is 
of the same order as the driving amplitude Because of these oscillations the tip-surface 
distance varies which then leads to a temporal modulation of the resonance frequency. Thus 
parametric resonance can take place. As an instability mechanism, parametric resonance is a 
threshold phenomenon. The cantilever starts to oscillate only if the modulation amplitude is 
larger than some threshold. Moreover, this threshold depends on the ratio of the modulation 
frequency and the eigen frequency. Its absolute minimum occurs when the modulation 
frequency is twice the resonance frequency, the so-called first-order parametric resonance 
condition. The minimal threshold is proportional to the damping constant (i.e., the inverse 
quahty factor) [Q. 

For a fixed modulation amplitude above the minimal threshold and a varying modula- 
tion frequency, one finds a frequency regime where parametric resonance occurs. As shown 
below, the center of the interval is determined by the second derivative of the interaction 
potential whereas the third derivative determines the width of the interval. If the mod- 
ulation amplitude is much higher than the minimal threshold, the effects of damping can 
be neglected. Then the maximum growth rate (which occurs roughly in the middle of the 
instability interval) is proportional to the third derivative of the potential. Because it does 
not depend on the Q factor, the instability may grow much faster than one would expect 
from the slow relaxation in a system with a high Q factor. 

After this qualitative discussion we calculate parametric resonance in an AFM by a model 
which describes the AFM by a single oscillator (see also Fig. |I|): 
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-^x ^ —x + Kx — V'{d + hcosujt — x)=Q, (1) 

where k is the spring constant, ujq = 27r/o is the resonance frequency, x is the tip position 
relative to its equihbrium value, d is the tip-surface distance in equilibrium without oscilla- 
tions, h and uj = 27r/ are the modulation amplitude and frequency, respectively, and V{r) 
is the tip-surface interaction potential. Derivatives of V are denoted by quotation marks, 
e.g., V'{r) = dV{r)/dr. 

The tip-surface interaction potential has in general an attractive and a repulsive part. 
In principle the proposed parametric resonance method can work in non-contact mode as 
well as in contact mode. Here we will restrict ourself to the non-contact mode (i.e., the 
tip is in the attractive part of the potential). In the intermittent or trapping mode ana- 
lytical approach becomes difficult because tip and sample deformations must be taken into 
account. Furthermore, the damping is more complicated because it increases significantly 
upon contact 0. 

In our model we use for the attractive part of the tip-surface interaction the following 
effective potential describing the sum of Van der Waals interactions between a plane 
sample and a tip with a spherical apex: 

= (2) 

where A is the so-called Hamaker constant, R is the tip radius, and r <S i? is the tip-surface 
distance. We have chosen A = 2.5 ■ lO^^^J and R = lOnm. 

In (P we have assumed that the tip-surface distance is modulated directly by vibrating 
the sample vertically. Alternatively one can excite the cantilever at its support which also 
leads to an oscillation of the tip-surface distance (see Fig. |l]). Here we restrict our discussion 
to the former case because in the single oscillator model both kind of driving are math- 
ematically equivalent. That is, driving the support of the cantilever with the amplitude 



H = hy[l — (cu/cjo)^]^ + {Q~^uj/u!o)'^ leads to a tip-surface oscillation with amplitude h. 

For h the equation of motion (|^) has a solution Xh(t) with the same periodicity as 
the modulation term, i.e., Xh{t + 2it/uj) = Xh{t) for any t. Parametric resonance means that 



the solution Xh becomes unstable. That is, a tiny perturbation Sx added to Xh increases 
exponentially in time. Linearizing (|1|) in 6x and assuming 



Xh<^h<^d (3) 



we get 



^5x + -^6x +{x + Vq + Vo'h cos ujt)6x = 0, (4) 



where 



and K'" = "^'^ 



(5) 

d 

This linearized equation of motion (|) is the well-known Mathieu equation with a damping 
term. The Floquet theorem implies that (^) has two linearly independent solutions ci{ujt)e^^^ 
and C2{ujt)e^^^ with ci/2(0 + 2it) = ci/2(0). The so-called Floquet exponent A is in general 
complex. Any solution of is a superposition of these solutions. The solution Xh is stable 
if any disturbance 6x decays. The stability condition is therefore ReAi/2 < 0. 

In order to calculate the Floquet exponents one has to expand the c functions into Fourier 
series. This expansion turns (^) into a set of infinitely many linear algebraic equations for 
the coefficients of the Fourier series. In the case of weak damping and weak driving, i.e., 

I V" I h V" 
Q>1, and ^^<1 + J^, (6) 

K K 

one gets good approximations for the Floquet exponents by restricting the Fourier series to 
the two leading components. That is, we solve (^) approximatively with the ansatz 

Sx = (a+e'"*/2 + a„e-'"*/2)e^* + c.c (7) 

Plugging the ansatz into and neglecting terms with exp{±3iujt/2) leads to two linear 
homogeneous equations for a± which have a nontrivial solution only if A is a solution of the 
characteristic polynomial 
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This equation can be solved analytically because it is a second order polynomial in (A/cJo + 
1/2Q)'^. Two solutions are always complex with a negative real part. The two other solutions 
become real if h is large enough. One of them is always negative whereas the other one 
becomes positives if 



h > hc = 
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where he is the threshold for parametric resonance. Typical growth rates given by the 
solution of d^) with the maximum real part is shown in Fig. |^. The threshold curve he and 
curves of constant growth rates are shown in Fig ^ 

The minimum value of h^ is given by 2Q~'^k/\Vq'\ which occurs at the parametric res- 



onance condition uj = 2ujq^1 + Vq / k. Note that the minimum value of he increase like 
d!^. Thus for increasing distance d the minimum value of the modulation amplitude that 
is necessary to lead to parametric resonance eventually exceeds the distance itself. Beyond 
that point which scales like Q^^^ parametric resonance is no longer possible. 

For h ^ 2Q^^k/\Vq"\ the damping term in (||) and (H) can be neglected. Using (|^) one 
finds that parametric resonance occurs for 



UJ^ < UJ < UJj, 



with u± = 2uJo\ 1 H ± 

V K 
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From measuring uj± one obtains Vq and Vq". The growth rate of a disturbance is zero at 



UJ = uj± and has its maximum near 2ojoJl + Vq / k. The maximum growth rate reads 



A. 
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Note that the handy formulas ( jlDD and (|TT]) hold only if the assumptions (0) and (^) are 
fulfilled. The parameters of Figs. ^ and |^ clearly fulfill these assumptions. 



In order to measure the (positive) growth rate in parametric resonance AFM we propose 
the following operation procedure. The procedure switches between two modes of operation: 
The on and the off mode. In the on mode the system is modulated at roughly twice the 
resonance frequency. In the response of the cantilever we are looking for oscillations at 
half of the driving frequency (i.e., near the resonance frequency of the cantilever) because 
such oscillations will be excited if parametric resonance occur. If there is an instability 
the amplitude a of these oscillations increases exponentially and eventually exceeds some 
predefined value 02. Then the system switches into the off mode where the modulation is 
switched off. The parametrically excited oscillation decays until the oscillation amplitude 
decreases below a value ai < 02 which is of the order of the noise level. Then the system 
switches back into the on mode. By measuring the average duration ton of the on cycles one 



gets the growth rate A. Assuming that x has to be increased from the noise level y (x^) up 
to 02 we get 

X ^ ^ '^^ . (12) 

^on 

Because the quality factor Q in vacuum is very high (i.e., of the order of 10^) the decay 
of the oscillation would take much more time than the excitation. In order to shorten toff 
one can use an external control circuit which brings the cantilever much faster back to 
equilibrium. 

We have simulated this operational scheme by integrating the equation of motion (|l]) 
numerically. In order to incorporate the effect of noise we have added on the right hand side 
of (|I]) an additive force 



with white noise, i.e.. 



(Kt))=0, {u{t)v{t + T))=5{r). (14) 

To extract the amplitude of the cantilever at half of the driving frequency we used a simple 
band filter, i.e., 



ij + 2-fuJoy + l^-j y = x. 




(15) 



For the relative bandwidth we have chosen 7 = 0.05. From the solutions of ( ]T5| ) we define 
an oscillation amplitude 



For x{t) = aocos{ujt/2) one gets exactly a = oq. In our simulation we have modeled 
the external control in the off mode by adding on the right-hand side of (|ID the force 
—Kx/{u!oQo)- Figure § shows an example of such a simulation. 

In Fig. H the growth rates obtained from such simulations are compared with the ana- 
lytical result, each square (with error bars) are the average of 30 to 50 on cycles. There 
is a good agreement between the analytical result and the simulation concerning uj±. The 
growth rate from the simulation is in general to small because Eq. ( p!2D assumes that the 



parametrically excited mode has initially the amplitude y (x^) neglecting the fact that in 
the noise there are also contributions from the stable modes. Thus the initial amplitude of 



With our simulations we have shown that parametric resonance AFM is possible. In- 
stability intervals and growth rates can be measured by using our proposed on-off scheme. 
From this data one gets with the help of ([l0| ) and ([TT|) the second and the third derivative 
of the tip-surface interaction potential. This method should work very well in the non- 
contact mode. Using this method in imaging we expect sharp contrasts because parametric 
resonance is a threshold phenomena. 
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the parametrically excited modes is on average less than w(x^). 
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FIGURES 
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FIG. 1. A schematic sketch of an atomic force microscope (AFM). The tip-surface distance is 
either modulated directly (modulation amplitude h) or indirectly by exciting an oscillation in the 
cantilever (driving amplitude H). 
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FIG. 2. Growth rate A as a function of the modulation frequency oj. Frequencies and growth 
rates are measured in units of the resonance frequency loq of the cantilever far away from the 
surface. The solid line is the solution of the characteristic polynomial (^). Squares show the result 
obtained from simulations as shown in Fig. ^. The parameters are d = 5nm, k = IN/m, Q = 10^, 
and h = 0.5nm. 
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FIG. 3. Instability threshold he and curves of constant growth rate A as functions of the 
modulation frequency uj. The parameters are the same as in Fig. 0. 
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FIG. 4. Simulation of an AFM. Solid lines show the on mode where the tip-surface distance is 
modulated with the frequency uj = 1.993u;o- This modulation is switched off when the amplitude 
a [defined by (p!6|)] reaches a threshold a2 (here: a2 = 0.6). In the off mode (dotted lines) the 
quality factor Q is set to a relatively low value Qq (here: Qq = 100). This simulates in a simplified 
way a control which brings the oscillating AFM back to equilibrium on a fast time scale. When 
the amplitude a is below a critical value ai < 02 (here: oi = 0.01) the system switches back to the 
on mode. The parameters are the same as in Fig. |. Noise level: (x^) = 10"^ 
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